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ABSTRACT 
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<^ ' We examine the reach of a Beta-beam experiment with two detectors at carefully chosen 
Q ; baselines for exploring neutrino mass parameters. Locating the source at CERN, the two detectors 
O and baselines are: (a) a 50 kton iron calorimeter (ICAL) at a baseline of around 7150 km which 
^ ' is roughly the magic baseline, e.g., ICAL@INO, and (b) a 50 kton Totally Active Scintillator 
^ . Detector at a distance of 730 km, e.g., at Gran Sasso. We choose '^B and '^Li source ions with a 
00 boost factor 7 of 650 for the magic baseline while for the closer detector we consider ^^Ne and ^He 
^ ^ ions with a range of Lorentz boosts. We find that the locations of the two detectors complement 
' each other leading to an exceptional high sensitivity. With 7 = 650 for '^B/^Li and 7 = 575 for 
^ ■ ^^Ne/^He and total luminosity corresponding to 5 x (1.1 x 10^^) and 5 x (2.9 x 10^^) useful ion 
decays in neutrino and antineutrino modes respectively, we find that the two-detector set-up can 
probe maximal CP violation and establish the neutrino mass ordering if sin^ 2^13 is 1.4 x 10~^ and 
2.7 X 10~^, respectively, or more. The sensitivity reach for sin^ 2^13 itself is 5.5 x lO"'^. With a 
factor of 10 higher luminosity, the corresponding sin^ 2^13 reach of this set-up would be 1.8 x 10~^, 
4.6 X 10"^ and 5.3 x 10~^ respectively for the above three performance indicators. CP violation can 
be discovered for 64% of the possible 6cp values for sin^ 26^13 > 10^^ (> 8 x 10"^), for the standard 
luminosity (10 times enhanced luminosity). Comparable physics performance can be achieved in 
a set-up where data from CERN to INO@ICAL is combined with that from CERN to the Boulby 
mine in United Kingdom, a baseline of 1050 km. 
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1 Introduction 



Long baseline experiments with powerful neutrino beams [1, 2, 3, 4] from upgraded accelera- 
tor facilities are the next frontier in neutrino oscillation physics. Set-ups with intense neutrino 
beam sources and smart detector technologies are being planned to fathom the hitherto uncharted 
regimes of the neutrino mass matrix. In particular, the long baseline experiments are being de- 
signed to measure the third mixing angle ^13, CP phase Sep and sgnlAml^) aka, the neutrino 
mass hierarchy-*- . The most promising avenue for this piirposc is the z/g oscillation channel 

Pg^j (or its T-conjugate -P^e), often referred to in the literature as the "golden channel" [5]. It was 
realized that while this channel can be used in the most cost-effective way to measure all three 
neutrino oscillation parameters listed above, it is also rife with the so-called problem of "parameter 
degeneracies". These are identified as the (^13, 5cp) intrinsic degeneracy [6], the (sg'n(Am|i), 5cp) 
degeneracy [7], and the {62^,1^/'^ — ^23) degeneracy [8]. Together they could result in as many as 
eight-fold degenerate solutions [9], of which, obviously, only one is true. This is evidently a very 
undesirable situation and a large body of existing literature is devoted to finding ways of combat- 
ing this menace [6, 10, 11, 12, 13, 14, 15, 16, 17]. A particularly attractive way of kiUing the clone 
solutions from the 5cp dependent degeneracies is to perform the experiment at the "magic base- 
line" [18, 19, 20]. At this baseline all 5cp dependent terms drop out, providing an ideal bedrock 
for measuring ^^13 and s(7n(Am|^). 

In [21, 22, 23] we expounded the sensitivity reach of a magic baseline experimental set-up 
where the neutrino source is a Beta-beam [24] at CERN and the detector is a large magnetized 
iron calorimeter (ICAL) at the India-based Neutrino Observatory^ (INO) [25]. Using a 50 kton 
fiducial mass and 80% detector efficiency^ for 1CAL@IN0, and considering a Beta-beam using 
*B and ®Li as ion source [26, 27], with a Lorentz boost 7 = 650 and assuming 5 x (1.1 x 10^^) 
and 5 x (2.9 x 10^^) useful ion decays in neutrino and antineutrino modes respectively, we showed 
that this set-up could unambiguously probe sgn{Am\i) at 3(T if*^ sin^ 26'i3 (true) > 5.6 x 10^^ and 
find a signal for ^13 at the same level of significance if sin^ 2^13 (true) > 5.1 x 10"^, independent 
of the true neutrino mass hierarchy and Sep (true). The ^13 and hierarchy sensitivity reach of this 
"magical" set-up is therefore superior to that of most other rival proposals involving Beta-beams 
[16, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39], and is almost comparable to that possible with 
a Neutrino Factory [40] at the magic baseline. This tremendous sensitivity of the CERN-INO 
Beta-beam project to ^13 and hierarchy comes from a combination of being close to the magic 
baseline as well as from being sensitive to near-maximal matter effects. 

The main drawback of the CERN-INO Beta-beam set-up is its inability to determine 6cp, since 
being almost at the magic baseline it is insensitive to the CP phase. Note that this is also true for 
the Neutrino Factory experiment at the magic baseline. This is why the optimal Neutrino Factory 
set-up demands a combination of two baselines, one magic and the other around L — 4000 km, 
where one has the best sensitivity to Sqp [40]. For Beta-beams, one has to optimize not only 

^In this paper -we define Amf^- = mf — and refer to sgn{Am'^i) as the neutrino mass hierarchy - sgn{Am'^i) > 
is called "normal hierarchy" (NH) -while sgn{Am'^i) < is called "inverted hierarchy" (IH). We stress that the 
arguments are valid for both hierarchical as -well as quasi-degenerate neutrino mass spectra. 

^The CERN to INO distance is 7152 km, which is tantalizingly close to the magic baseline. 

^We give details about the detector configuration in section 3. 

^Thcsc are the values of the parameters chosen by Nature, to be distinguished from the fitted values. Throughout 
this paper we denote the true value of a parameter by putting "(true)" after the symbol for the parameter. 
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on 7, the baseline L, and the luminosity, but also on the choice of the source ions. A detailed 
optimization study has revealed that for intermediate values of 7, one would preferentially use 
and ^Li at the magic baseline for the mass hierarchy, ^^Ne and ^He at an intermediate basehne 
of L = 600 — 900 km for CP violation, and either ^^Ne and ^He at the intermediate baseline or 

and ^Li at the magic baseline for the sin^2^i3 discovery [41]. 

In this paper we consider a two-baseline Beta-beam set-up, one with L = 7152 km, the CERN- 
INO basehne, and another with L — 730 km which is the CERN-Gran Sasso (LNGS) distance. 
For the CERN-INO case ^B and ^Li are the preferred source ions and we take 7 = 650. For the 
CERN-LNGS set-up, on the other hand, we choose the ^^Ne and ®He ions and allow their 7 to 
vary between 250-650^ in order to examine the dependence of the sensitivity on the value of the 
boost factor. Since the ^B and ®Li ions would produce multi-GeV neutrino beams for 7 = 650, 
we use a 50 kton iron calorimeter (ICAL) for the longer basehne at INO and call this set-up 
CERN-ICAL@INO. For the intermediate basehne option, since we are interested in the lower 
energy ^^Nc and ^He ions, we assume a 50 kton Totally Active Scintillator Detector (TASD) in 
order to harness the low energy events required for better CP sensitivity. We present results for 
5 X (1.1 X 10^^) and 5 x (2.9 x 10^*) useful ion decays in neutrino and antineutrino modes respectively 
for both baseline set-ups. We also show the projected sensitivity if one achieves statistics that are 
larger by a factor of 10. Since the statistics depend on the product of the size of the detector, the 
exposure time, detector efficiency and the number of useful ion decays, this one order of magnitude 
increase could come from a combination of enhanced performance in any of the factors mentioned 
above. In particular, it might be possible to have 10 times larger useful ion decays per year in the 
storage ring [43]. ICAL@INO could also be upgraded to 100 kton. 

Note that while most of our figures and discussion in the text would explicitly be for the 
CERN-LNGS intermediate baseline set-up, we also present results for the case where a ^^Ne and 
^He Beta-beam is shot from CERN to the Boulby mine in the United Kingdom (see [38] and 
references therein). The CERN-Boulby distance is 1050 km, and one can check from Fig. 5 of 
[41] that the CP sensitivity for this baseline is only marginally weaker than that for L — 730 
km. Therefore, when combined with the CERN-INO data, we expect similar performance for this 
baseline option as well. 

In some sense this work is a part of an ongoing international exercise in sharpening the full 
capability of long baseline Beta-beam experiments to explore neutrino properties. In addition 
to the work discussed earlier in the Introduction, some of the milestones on this route can be 
identified as the CERN-MEMPHYS proposal [16, 28, 29], its variants with higher 7 [32, 33], using 
a cocktail of sources [37], detector options, etc. 

While the best CP violation sensitivity comes with ^^Ne and ^Hc as source at our chosen 
intermediate baseline, it might be more convenient and less demanding to use the same set of 
ions at both the basehnes. The optimal basehne for CP studies with ^B and ^Li was seen to be 
around 1000-2000 km in [41]. In [31], the authors proposed a set-up with this scenario where they 

^With the existing facilities at CERN, one will be able to accelarate ^He to only about 7 = 150. This corresponds 
to a boost factor of about 250 for ^*Ne and 280 for *B, due to the different charge to mass ratios for these ions. 
However, with the "Super-SPS", an upgraded version of the SPS with super-conducting magnets [33, 30] it should 
be possible to accelarate ®He to 7 = 350, which corresponds to 7 = 575 for -'^^Ne and 7 = 650 for *B. The Tevatron 
at Fermilab could produce a Beta-beam with similar boost factors. Higher acceleration at CERN would require 
the use of the LHC itself, with 7 > 1000 possible [42]. 
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considered a two-baseline combination with a longer L ~ 7000 km and a shorter L ~ 2000 km. A 
7 = 350 Beta-beam produced by and ^Li as source ions feeds both detectors and they consider 
up to 5 X (10 X 10^^) total useful ion decays in both neutrino and antineutrino modes. They 
consider 50 kton of magnetized iron as the end detector for both baselines. In this set-up the mass 
hierarchy can be discovered if sin 26*13 (true) > 3 x 10"^ and CP violation can be established for 
70% of the possible 5cp (true) values if sin^ 20i^ (true) > 10~^. Signal for the mixing angle itself 
could be observed for sin^ 2^13 (true) > 1 x 10~^. For 5 x (2 x 10^^) total useful decays, which is 
comparable to our standard luminosity, the corresponding reaches are sin^ 2^13 (true) ^ 10~^ for 
the mass hierarchy, and sin^ 2^13 (true) ^ 8 x 10~^ for discovering ^13. CP violation sensitivity 
gets severely restricted and deteriorates, especially for 5cp (true) = 270°. The poor CP reach 
for lower luminosity comes due to the choice of ^B and ®Li as source ions (see [41] for a detailed 
discussion) . Since much better CP sensitivity for plausible luminosity and 7 can be achieved with 
the ^^Ne and ^He combination, we will exclusively use them for our intermediate baseline option 
of CERN to LNGS. Since our chosen 7 is larger, wc also achieve sensitivity to mass hierarchy and 
^13 which is better by at least an order of magnitude for the same number of useful ions decays. 

The paper is organized as follows. In the next section we briefly discuss the oscillation proba- 
bility for the golden channel and highlight the magic baseline behavior. The experimental set-up 
which we consider is outlined in the subsequent section. In the following section we explore the 
prospects of the TASD@LNGS detector, in a stand-alone mode, for exploring CP violation. This 
is to set the stage for the next section where the results of a combined two-baseline analysis are 
presented with emphasis on the sensitivities to CP violation, 6*13, and the neutrino mass hierar- 
chy. We end by summarizing the results for our combined CERN-INO and CERN-LNGS set-ups 
and make a comparative discussion of this set-up vis-a-vis the combined CERN-INO and CERN- 
Boulby set-up. We also compare our set-up against the one studied in [31] and discuss the pros 
and cons of the two proposals. Some remarks about the beam related as well as atmospheric 
neutrino backgrounds are collected in the Appendix. 



For the results presented in this paper we have calculated the neutrino oscillation probability in 

matter exactly, using the PREM profile for the Earth matter density [44] . However, for elucidating 
the behavior of neutrino oscillations as a function of baseline and/or neutrino energy, it is useful 
to exploit the approximate analytic formula for the golden channel probability Pefi. in matter [45, 
46, 47], keeping terms only up to second order in the small quantities 6*13 and a = Amli/ Am'^^ 



2 Golden channel oscillations 



[5, 48] 
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4 



+ a'cosH2ssm'2eu ^'{'^^\ (1) 



where 



Am^i L ^ A 

^31 



^=±2^G.A'.£. (2) 



Above, A is the matter potential, expressed in terms of the electron density A^e and the (anti) neutrino 
energy E; the '+' sign refers to neutrinos while the ' — 'to antineutrinos. 
Prom Eq. (1) it is apparent that at a baseline for which 

sin(iA) ~ , (3) 

the oscillation probability is insensitive to 6cp- This defines the "magic baseline". As noted, the 
CERN-INO distance closely matches this baseline and it has been shown in earlier work [21, 22, 23] 
that this ensures a high sensitivity at the same time to the determination of 9is and the extraction 
of the neutrino mass hierarchy. In this work, we examine the benefits which accrue from combining 
two Beta-beam experiments both with the source at CERN, one with the detector at INO and the 
other located at an appropriate distance where the impact of CP violation is prominent. Taking a 
cue from an optimization analysis of such Beta-beam set-ups [41], this second detector is assumed 
to be at Gran Sasso. 

We will refer to Eqs. (1) and (2) from time to time to identify the physics underlying the 
results which we present later. 



3 The Experimental Set-up 

The number of (anti)muon events in the ith energy bin in the detector is given by 

Tn„ fine /-Smax /-Br" 
Ni ^ ; ;y / dE : dEAHE)a.,{E)R{E,E^)P,,{E), (4) 
47rL^ Jo JET'' 

where T is the total running time, n„ is the number of target nucleons in the detector, fm is the 
charge identification efficiency (meaningful only for the magnetized iron detector) , e is the detector 
efficiency and R{E, Ea) is the energy resolution function of the detector, which we assume is a 
Gaussian. Por muon (antimuon) events, a^^ is the neutrino (antineutrino) interaction cross- sect ion. 
The quantities E and Ea are the true and reconstructed (anti)neutrino energies respectively. The 
quantity (t){E) is the (anti)neutrino Beta-beam fiux produced at the source, with Lorentz boost 7 
and L is the baseline. 



3.1 The Source Ions 

The technological challenge for producing the Beta-beam involves creating, bunching, accelerating 
and storing beta unstable radioactive ions [49, 50]. A pure and intense v^. and/or Pg beam is 
produced when these highly accelerated ions decay along the straight sections of the storage ring. 
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Ion 


r(s) 


£^0 (MeV) 


/ 


Decay fraction 


Beam 




2.41 


3.92 


820.37 


92.1% 






1.17 


4.02 


934.53 


100% 




IB 


1.11 


14.43 


600872.07 


100% 




lu 


1.20 


13.47 


425355.16 


100% 





Table 1: Beta decay parameters: lifetime r, electron total end-point energy E'o, /-value and decay 
fraction for various ions [51]. 



The choice of ions is governed by a number of factors which include their end-point energy, life- 
time, production rate and their charge to mass ratios. Two sets of source ions have been identified 
as plausible candidates for the production of a Beta-beam. The most widely discussed ions are 
^^Ne and ^He for z/g and z/g beam respectively. The alternative set of ions which have also been 
considered extensively in the recent literature are ^B and ^Li for u^. and u^. beam respectively. We 
present in Table 1 the relevant details of the properties of these ions. We note that the main 
difference between the two sets lies in their total end-point energies - ^B (^Li) has an end-point 
energy, Eq, which is 3.68 (3.35) times that of ^^Ne (^He). The Beta-beam flux spectrum depends 
solely on 7 and Eq and the energy of the beam is often crucial in determining the type of physics 
that can be explored. We therefore reiterate the following points [41, 4]: 

• Since the peak neutrino energy is roughly given as 7-E'o) the boost factor needed for ^^Ne (^He) 
should be about 3.68 (3.35) times that needed by ^B (^Li) in order to achieve the same 
neutrino peak energy. 

• Since the total number of neutrinos peaked in the forward direction, i.e., towards the de- 
tector, roughly goes as 7^, larger 7 enhances the flux. Higher boost factors are therefore 
preferred. 

• The reference number of useful decays which we use for antineutrinos (2.9 x 10^^ per year) 
is larger than that for neutrinos (1.1 x 10^*^ per year). The ratio between the two just about 
compensates for the smaller antineutrino interaction cross section in matter. So, roughly, 
this results in a symmetry between the v and P data. 

We have seen in [22, 23] that the higher end-point energy of ^B and ^Li allows these ions to 
produce Beta-beams with peak energy in the multi-GeV regime, where one obtains near-resonant 

matter effects for the near-magic baseline. This can be achieved with plausible values of 7, 
and results in remarkable sensitivity to ^^13 and sgn{^m\^). Indeed this set of ions emerged 
as the preferred choice for the magic baseline after a thorough scan of plausible 7 values [41]. As 
mentioned before, the 5cy> sensitivity is smothered at the magic baseline. It was seen in [41] that 
intermediate baselines are far superior for CP studies. It was further noted that CP sensitivity 
with ^^Ne and ^He is much better than that with *^B and ^Li. Therefore, we take ^B and ^Li as 
source 10ns for the CERN-INO baseline and ^^Ne and ^He as the source ions for the CERN-LNGS 
baseline. We will fix the Lorentz boost of ^B and *^Li for CERN-INO as 7 = 650 and will study 
the impact of 7 for ^^Ne and ''He at the shorter baseline. 
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Detector Characteristics 


ICAL@INO 
(Only 


TASD@LNGS 
(Both ij,^ & e^) 


Total Mass 


50 kton 


50 kton 


Energy threshold 


1 GeV 


0.5 GeV 


Detection Efficiency (e) 


80% 


80% (/x±) & 20% (e±) 


Charge Identification Efficiency (fio) 


95% 


No CID 


Energy Resolution (cr) (GeV) 


0.15E(GeV) 


0.03^E(GcV) for 
0.06^E(GeV) for e± 


Bin Size 


1 GcV 


0.2 GcV 


Background Rejection 


0.0001 


0.001 


Signal error 


2.5% 


2.5% 


Background error 


5% 


5% 



Table 2: Detector characteristics used in the simulations. The bin size is kept fixed, while the 
number of bins is varied according to the maximum energy. 



Throughout this work we deal with results obtained by combining data from neutrino and 
antineutrino beams. Further, we ascribe the same Lorentz boost 7 to both beams. It has to be 
borne in mind that the charge to mass ratios, {Z/A), of the source ions are not the same (see 
Table 1), and so the same boost cannot be achieved if the ions run simultaneously in a ring. Note, 
however. 



Consequently, if the two neutrino (antineutrino) source ions are run simultaneously in the ring 
in which the ^B (^Li ) ions have a Lorentz boost 7 = 650, then the ^^Ne (^He ) source ions will 
have a 7 of 575. Of course, these are but sample boost factors for illustrating the feature that 
simultaneous run configurations are possible for the chosen source ions, with the proviso that the 
(^B and ^Li ) pair have a different value of 7 related to that of the (^^Ne and ^He ) pair. 



3.2 The ICAL@INO detector 

The ICAL detector at INO will be a 50 kton magnetized iron calorimeter with Resistive Plate 
Chambers (RPCs) serving as the active elements [25]. The first phase of this detector (starting 
around 2012) will focus on atmospheric neutrino measurements. The entire ICAL detector will be 
split into 3 modules, each of which would be of dimension 16 m x 16 m x 12 m. The detector 
structure consists of iron layers of 6 cm thickness with a 2 cm gap between layers wherein the 
glass RPCs will be interleaved. The detector would be magnetized by a field of about 1 Tesla. 
The detector characteristics that we have used for our simulation are shown in Table 2. It is 
planned that ICAL@INO will be upgraded to 100 kton in the future. Number of expected muon 
(or antimuon) events can be calculated using Eq. (4). 
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Figure 1: Total event rates expected in 5 years at the CERN-TASD@LNGS set-up assuming 
1.1 X 10^^ and 2.9 x 10^® useful decays per year for ^^Ne and ^He , respectively. The left panel 
shows the dependence on the boost factor at sin^ 26*13 = 10~^, while the right panel depicts the 
variation with sin^ 26*13 for a fixed 7 = 500. Results are shown for four values of 6cp- 0°, 90°, 180°, 
and 270° for both neutrinos (thick broken lines) and antineutrinos (thin broken lines). Normal 
hierarchy has been assumed. In the left panel, the vertical line corresponds to 7 = 575 (see text) 
and the solid curves are the effective neutral current backgrounds for neutrinos (thick red lines) 
and antineutrinos (thin red lines). 



3.3 The TASD@LNGS detector 

For the CERN-LNGS baseline we will be working with a 50 kton TASD. A scintillator detector has 
several virtues. It has good detection efficiency for muons, can support a low energy threshold and 
has excellent energy resolution. It can also detect electrons but with a lower efficiency. On the flip 
side, unlike in the iron calorimeter, there is no charge identiflcation. The detector characteristics 
we have assumed [2, 35] are listed in Table 2. One can see from the Table that the background 
suppression capability of a TASD is also expected to be one order poorer than that of an iron 
calorimeter detector^. The effect of this on the results is discussed later. However, the biggest 
advantage of TASD over iron calorimeters is its much lower threshold energy. This helps in 
eliminating the clone solutions and hence delivers better CP sensitivity. 

We begin with some remarks about the performance of the CERN-TASD@LNGS Beta-beam 
set-up in isolation before turning to our main theme, namely, a combined analysis of the results 
from two detectors at different baselines. The number of useful ion decays that we choose as 
reference values are (1.1 x 10^^) per year for ^^Ne and (2.9 x 10^®) per year for ^He . For comparison 
we also present results obtained using luminosities one order higher than these reference values. 
We combine data from u and u runs, and assume that both are at the same Lorentz boost 7. The 

^This can be traced to the fact that in scintillator detectors hadrons travel farther before they stop or decay, 
enhancing the possibility that they mimic the signal. The lack of charge identification is also a handicap in this 
regard. 
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values of the oscillation parameters and the earth's matter density profile that we use are listed 
in Table 3. 

In Fig. 1 the five-year (anti)muon event rates for the TASD detector at Gran Sasso are shown. 
In the left panel the number of events is shown as a function of 7 for sin^ 2^13 = 10~^ while in 
the right panel it is shown as a function of sin^ 2^13 when 7 is fixed at 575. In both panels results 
are shown for both v and v beams. Normal neutrino mass hierarchy has been assumed. Notice 
that in both panels, to a good approximation, the number of events for neutrinos for any 5cp 
closely matches that for antineutrinos but with CP phase —5cp- This behavior can be readily 
understood from Eq. (1). At the CERN-LNGS shorter baseline, matter effects are small and 
hence A « 1 in Eq. (2). It follows that the 5cp dependence of the probabilities can be expressed 
as Pen ~ ^0 + ^-sin^cp + A+cos^cp and Pg/i — ^0 — /l_sin5cp + /1+cos5cp, where Aq,A± 
are independent of 5cp, whence the symmetry is manifest. In the left panel, the increase in the 
number of events with increasing Lorentz boost 7 is a consequence of the resultant enhancement 
of the on- axis flux as well as the higher energy of the beam. The reason why neutrino event rates 
are higher for 5cp = 90° and 0° compared to that for 5cp = 180° and 270° can be easily seen 
from the above expressions. One can explicitly check that and /!_ are positive, which means 
that for 5cp — 180° and 270° number of events are given by a difference of two positive terms 
while for 5cp — 90° and 0° they come as the corresponding sum. The reason why the slope of the 
curves for 5cy> = 180° and 5cp = 270° (90°) for neutrinos (antineutrinos) is relatively smaller than 
that of 5cp = 0° and 5cp = 90° (270°) for neutrinos (antineutrinos) can be understood from a 
similar reasoning. The behavior of the number of events with sin^ 2^13 (right panel) is not difficult 
to explain. As sin^ 2^13 increases from a near-zero value, first the 2nd and 3rd terms in Eq. (1), 
which are linear in sin 2^13 and are dependent on 5cp, begin to contribute. Beyond a certain point, 
the first term takes over leading to the increasing behavior for all curves in the right panel. 

For the purpose of illustration, in the left panel of Fig. 1 we also present the v and v NG 
background estimates. Notice that the background is quite significant for the entire range of 7. 
We have also indicated the point corresponding to 7 = 575, which, as discussed earher, is the 7 
for ^^Ne (^He) ions when they are run in the same storage ring as (^Li) with the 7 for the latter 
being 650. 

Note that in our analysis for this detector, we have also included the information from the 
Pee and Pgg channels, albeit with the reduced detector efficiency (sec Table 2) of only 20% for 
electrons. The number of electron events can be calculated using Eq. (4), by making appropriate 
changes to the oscillation probability and cross-sections. However, as was noted in [13, 14], this 
channel has hardly any sensitivity to ^13 and mass ordering at low baselines like the one under 
discussion. It is independent of 5cp- 

4 CP sensitivity: CERN-TASD@LNGS Set-up 

The goal of this work is to explore the power of a two-detector Beta-beam experiment to unravel 
neutrino mass parameters. The detector at the magic baseline is insensitive to the GP phase 5cp 
and so the main emphasis of the closer detector will be to address this issue. To set the stage for 
the full two-detector analysis, in this section we discuss the GERN-TASD@LNGS combination as 
a stand-alone set-up for exploring GP violation. 
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Benchmark Values 


la estimated error 


Am^i(true)| = 2.5 x 10"^ eV^ 


a(|Amii|) = 1.5% 


sin^ 2^23 (true) = 1.0 


cT(sin" 2^23) = 1% 


Am^i(true) = 8.0 x 10"^ eV^ 


a(Amii) = 2% 


sm^^i2(true) = 0.31 


a(sm" Ou) = 6% 


p(true) = 1 (PREM) 


a{p) = 5% 



Table 3: Chosen benchmark values of the oscillation parameters and the earth matter density 
profile and their la estimated error. 



The 3a sensitivity to "maximal CP violation" is presented in Fig. 2 as a function of the boost 
factor 7: these are the (red) solid curves marked "CG" in both panels. This performance indicator 
is defined as follows. We generate the data for 5cp (true) = 90° (left panel) or 6cp (true) — 
270° (right panel), while in theory we allow Sep to be 0° and 180°. We marginahze over the 
normahzation of the PREM density profile and all other oscillation parameters, including the 
neutrino mass hierarchy. The lowest value of sin^ 2^13 (true) which allows us to rule out at 3a a 
CP conserving theory (both 0° and 180°) when CP is maximally violated in Nature, is plotted in 
the figure as the sensitivity reach for maximal CP violation. 

The close similarity of the sensitivity reaches displayed in the two panels of Fig. 2 is in line 
with expectation^. Prom the approximate symmetry 

-Pe/i < — ^ Pep under Sep < — > -Sep (6) 

noted earlier from Eq. (1), and recalling that the data involves u as well as P beams, it is to be 
expected that results for Sep (true) = 90° and Sep (true) = 270° will be almost identical^. More 
explicitly, the statistics for Sep (true) = 270° is comparatively poor for neutrinos with NH true, 
which is made up by the higher statistics for antineutrinos. The effect of the relatively lower 
statistics for neutrinos shows up mainly as additional degeneracies. As discussed in the previous 
section, the better energy resolution and lower threshold of the TASD results in an amelioration 
of the correlation and degeneracies and the islands appearing in Fig. 5 of [41] are absent here. 

We have assumed NH to be the true mass hierarchy in these plots. This is not expected to be a 
serious issue as at this baseline matter effects are small. In fact, the symmetry in Eq. (6) becomes 
exact if in addition to Sep ^ —Sep the neutrino mass hierarchy is also flipped. In view of this, for 
IH true, the left (right) panel of Fig. 2 would match with Sep (true) = 270° {Sep (true) = 90°). 
Of course, this argument rests on Eq. (1), which is an approximation, while the results are based 
on the exact expressions. We have checked that indeed taking IH to be true gives very similar 
results. 

It is seen from Fig. 2 that in both panels the sensitivity reach for maximal CP violation 
improves sharply with 7 initially, and then it more or less saturates beyond 7 ^ 500. This can be 

^This feature was also noted in [41] where the detector was a 50 kton magnetized iron calorimeter. 
^This symmetrical behavior breaks down for larger values of sin^ 2^13 (true). See later. 
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Figure 2: The 3a sin^ 2^13 (true) reach for sensitivity to "maximal CP violation" . The red 
solid curves (marked as 'CG') are for the CERN-TASD@LNGS alone while the blue dashed 
lines (marked as 'CG+CI') are for the combined data from CERN-TASD@LNGS and CERN- 
ICAL@INO. The results are shown as a function of the Lorentz boost for ^^Ne and ^He (taken 
same for both ions), for Sep (true) = 90° (left panel) and Sqp (true) = 270° (right panel). Thick 
lines (marked "High") are for 5 x (1.1 x 10^^) useful ^^Ne and ®B decays and 5 x (2.9 x 10^^) useful 
^He and ^Li decays, while thin lines (marked "St") are for 5 x (1.1 x 10^^) and 5 x (2.9 x 10^^) 
useful ion decays respectively. 



understood as follows. For L = 730 km, P^^ peaks around 1.5 GeV, and is very small beyond ~ 3 
GeV. For ^^Ne and ^He ions, the peak (anti) neutrino energy reaches 1.5 GeV for 7 ~ 430. For 
7 > 430, while the peak energy shifts linearly to higher values and the total flux increases as 7^, 
the flux at the energy ii^ ~ 1.5 where Pen is high remains roughly the same (cf. Fig. 1 of [41]). 
Therefore, for all 7 > 430 we do not see any marked improvement. 

In Fig. 2 results are shown for the standard luminosity as well as for the case where the 
total statistics is increased by a factor of ten. Increasing the luminosity has a significant impact 
on the sensitivity of the experiment. From the figure one can read off that for 7 = 575 and 
Scp (true) = 90°, CP violation can be discovered at 3cr for sin^ 2^13 (true) > 1.6 x 10~^ for the 
standard luminosity and for sin^ 2^13 (true) > 1.97 x 10""^ if luminosity is increased by a factor of 
ten. For (5cp (true) = 270°, these limits are only slightly worse, namely, sin^ 2^13 (true) > 1.8x10"^ 
and sin^ 26'i3 (true) > 2.03 x 10~^, respectively. Note that the results are for the NH true case 
but, as discussed earlier, for IH true the numbers will hardly change. 

In Fig. 3 we show the impact of changing the detector characteristics on the maximal CP 
violation sensitivity of the TASD@LNGS set-up. The sensitivity is shown in this figure as a 
function of 7 for a 5-year run with the reference choices of the {u, v) luminosities. The left panel 
shows the effect of changing the background fraction from the standard assumed value of 10~^ (red 
solid lines). We have shown for comparison results for background fractions of 10~^ (black dashed 
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Figure 3: Effect of changing the detector characteristics on the "maximal CP violation" sensitivity. 
The left panel shows the impact of changing the background fraction to one-tenth or ten times 
the reference value of 10~^ (red solid curves). Results are shown for 6cp (true) = 90° by thick 
lines and Sqp (true) = 270° by thin lines. The location of 7 = 575 is also shown. The right 
panel shows the effect on changing the detector threshold. Results for 6cp (true) = 90° are shown 
in the inset. The red solid curve corresponds to a threshold of 500 MeV while the blue dashed, 
black-dot-dashed, and pink triple-dot-dashed curves correspond to thresholds of 1 GeV, 1.25 GeV, 
and 1.5 GeV, respectively. 



lines) and 10~^ (blue triple-dotted-dashed lines). We conclude (at 7 = 575) that while lowering the 
background fraction by a factor of 10 improves the sensitivity by a factor of 1.7, a 10-fold increase in 
the background would deteriorate the sensitivity by a factor of 6. Right panel of the figure shows 
the impact of changing the detector threshold on the CP sensitivity. Increasing the threshold 
could help us in reducing the backgrounds^. We see that for 6cp (true) = 270° increasing the 
detector threshold results in the appearance of degenerate solutions. For 6cp (true) = 90° (inset) 
the effect is less severe; however, the sensitivity falls for both choices of 6cp (true) as the threshold 
is increased. The appearance of clone solutions (multiple solutions for the same 7 at different 
values of sin^ 2^13 (true)) when the threshold is increased beyond 1 GeV reflects the important 
role of lower energy [i.e., lower oscillation wavelength) neutrinos in this analysis. A strength of 
the scintillator detector is the scope of inclusion of the low energy data by virtue of the permissible 
lower threshold. We have also studied the effect of changing the bin size of the data on the CP 
sensitivity of our experimental set-up. We find the bin size as large as 1 GeV could also be 
accommodated without significantly deteriorating the CP violation sensitivity of the experiment. 

Fig. 2 gives the "maximal CP violation" sensitivity reach, i.e., the ability to distinguish 
6cp (true) = 90° or 270° from no CP violation using the CERN-TASD@LNGS detector. These 

^The atmospheric neutrino background would be larger at lower energies. Even backgrounds from neutral 
current events could be reduced by going to higher threshold. 
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Figure 4: The area enclosed by the curves represents the 3a range of 6cp (true) as a function of 
sin^ 2^13 (true) for which the data can be used to rule out the CP-conserving scenario using the 
CERN-Gran Sasso reference TASD set-up with ^^Ne and ^He as source ions. The Lorentz boost 
(the same for both ions) is fixed at 575. Results are shown for 5 years run with the reference 
luminosity (blue dashed curves) as well as for a luminosity which is one order higher (red solid 
curves) . 



choices of 6cp (true) are the extremal ones. What about other values of 6cp (true)? In Fig. 4 we 
show for 7 = 575, the range of 6cp (true) as a function of sin^ 2613 (true) for which one can use 
the data to rule out no CP violation at the 3a level. We show results with 5 x (1.1 x 10^^) and 
5 X (2.9 X 10^^) useful ion decays in neutrino and antineutrino modes respectively (the blue dashed 
curves) as well as those with one order of magnitude enhanced number of events in both channels 
(red solid curves). At the sensitivity reach values of sin^2^i3 (true) displayed in Fig. 2, this range 
is reduced to two points. The range grows rapidly as sin^ 26'i3 (true) increases to about 1 x 10^'^ 
beyond which it saturates, except for the kinks in the third and fourth quadrants of 6cp (true) 
around sin 2^13 (true) = 1x10-2, which arise from clone solutions. For sin 2^13 (true) > 10 ^ we 
should be able to discover CP violation for 64% of the possible 6cp (true) values for the standard 
luminosity. If luminosity was enhanced by a factor of 10, this would improve to 88% of the possible 
6cp (true) values. In fact, for 10 times larger statistics, CP violation can be discovered for 64% of 
the possible ^cp (true) values for Oi^ as low as sin^ 26*13 (true) > 8 x 10"^. 
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Figure 5: The left and right panels depict the sgnlAm^i) sensitivity reach and the sin^2^i3 
sensitivity reach, respectively, at 3a as a function of the boost factor for ^^Ne and ^He. In both 
panels, the red solid lines are for CERN-TASD@LNGS alone while the blue dashed lines are for the 
combined data from CERN-ICAL@INO and CERN-TASD@LNGS. Resuhs for 6cp (true) = 7r/2 
and 37r/2 are shown. Thick lines are for 5x (1.1 x 10^^) useful ^^Ne and decays and 5x (2.9 x 10^^) 
useful ^He and ^Li decays, while thin lines are for 5 x (1.1 x 10^*) and 5 x (2.9 x 10^^) useful ion 
decays respectively. The sensitivity reaches for the CERN-ICAL@INO set-up alone are indicated 
for both luminosities (for 6cp (true) = in the left panel) by arrows on the right side of the panels. 
The location of 7 = 575 is shown. 



5 Two-baseline combined results 

We are now in a position to discuss the benefits accruing from combining the results of the two 
detectors: one at a near-magic baseline distance and the other located 730 km from the source. We 
discuss in turn the sensitivity reaches of this combined set-up to maximal CP violation, sin^26'i3, 
and the mass hierarchy. 

Recall that the very long baseline detector is a magnetized iron calorimeter at the India-based 
Neutrino Observatory at a distance of 7152 km from CERN. For this detector we consider and 
^Li as the source ions for u and z/ beams and we assume that both are accelerated to a reference 
Lorentz boost 7 = 650. The reference choice for the number of useful *B (®Li ) ion decays is 
5 X (1.1 X 10^*^) (5 X (2.9 X 10^^)). We also show results for a 10 times enhanced luminosity. 

5.1 CP Sensitivity of the Combined Set-up 

We first probe the impact of adding the data from the CERN-ICAL@INO set-up on the sensitivity 
reach in sin^ 2^13 (true) for maximal CP violation. We chose to use the CERN-LNGS baseline with 
^^Ne and ^He as source ions since it was shown in [41] to be the ideal set-up for probing 6cp- In Fig. 
2 is shown (blue dashed curves) the impact of adding the CERN-ICAL@INO ®B and ®Li results 
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on the sensitivity to maximal CP violation. It is seen that there is insignificant improvement in 
the sensitivity. We can conclude that the CERN-LNGS basehne, which is the optimal choice as a 
stand-alone set-up for CP measurements, remains so after the addition of the new magic baseline 
data. This should not be surprising since at the magic baseline there is no Sep dependence, and the 
CP violation measurement of the 730 km experiment is almost devoid of parameter degeneracies. 

5.2 The sgn(Am|j^) Sensitivity Reach 

Next, we consider the potential of the two-basclinc Beta-beam set-up to explore the nature of the 
neutrino mass hierarchy. The "s5'n(Am|^) sensitivity reach" is defined as the range of values of 
sin^ 26*13 (true) for which the experiment could rule out the wrong hierarchy at 3a. We show this 
as a function of the boost of the ^^Ne and ^He ions in the left panel of Fig. 5. The boost for 
the *B and ^Li source ions for CERN-ICAL@INO is fixed at 7 = 650. We have shown results 
for 6cp (true) = 90° and 270° and assuming NH to be true. The solid lines show the sgn{Am^^) 
sensitivity of CERN-TASD@LNGS alone while the dashed lines give the sensitivity when the two 
data sets are combined. The sgn{Am'^i) sensitivity reaches for the CERN-1CAL@1N0 set-up with 
7 = 650 with the standard and enhanced luminosities (for 5cp (true) = 0) are indicated in the 
figure by arrows. In Table 4 we give for a fixed 7, the sensitivity reach in sin^ 2^13 (true) for 
determining the mass hierarchy for the single and combined baseline set-ups considered in this 
paper. A comparison with the sensitivity reach envisaged for the optimized two-baseline neutrino 
factory from Ref. [4] is also given for comparison. From the figure and the table we note that the 
hierarchy sensitivity is rather poor for the CERN-TASD@LNGS set-up alone since the baseline is 
not large enough to capture enough matter effect and is almost independent of the value of 7. As 
shown in [23], the CERN-ICAL@INO experiment has a very good hierarchy sensitivity. Notice the 
near-independence on 6cp (true) - a feature of the magic baseline. For the combined two-baseline 
set-up we achieve an exceptional sensitivity such that (for 6cp (true) = 90°) the mass hierarchy 
could be determined at 3a if sin 2^13 (true) > 2.7 x lO"'' and 4.64 x 10"^ for the standard and 
ten times enhanced luminosity respectively. All results presented are for NH true, however, we 
have checked that results with IH true are similar. 

5.3 The sin^ 2^13 Sensitivity Reach 

Finally, we turn to the expected sensitivity of the combined set-up to ^13. We define the perfor- 
mance indicator "sin^ 2^13 sensitivity reach" as follows. We generate the data at sin^ 2^^i3 (true) = 
and determine the values of sin^ 2^13 that could fit this data within a chosen confidence level, 
allowing all oscillation parameters, including 6cp and the mass hierarchy, to vary freely in the fit. 
We also marginalize over the normalization of the PREM density profile. This is applicable when 
the experiment has observed no ^13 driven signal and at best can put an upper limit on the still 
allowed sin^ 2^13. The 3a projected sin^ 2^13 sensitivity as a function of 7 for ^'^Ne and ®He is shown 
in the right panel of Fig. 5. The boost of and ^Li for the CERN-ICAL@INO data is fixed at 
7 = 650. The red solid lines in the figure show the sin^ 2^13 sensitivity for CERN-TASD@LNGS 
alone, while the dashed lines are what we expect by combining the two data sets. The sensitivity 
reach in sin^ 26*13 for the CERN-ICAL@INO is given in Table 4 and is indicated in the figure by 
an arrow. As can be seen from Table 4, the sensitivity reach for CERN-TASD@LNGS alone is 
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comparable to that of CERN-ICAL@INO taken alone. For ten times enhanced statistics the cor- 
responding sin^ 2^13 sensitivity reach of CERN-ICAL@INO and CERN-TASD@LNGS taken alone 
are, sin^ 2^13 < 1.76 x 10~^ and sin^ 2^13 < 8.59 x 10~^ respectively. A combination of the two 
data sets improves the performance tremendously. For the standard luminosity we could reach 
down to sin^ 2^13 < 5.46 x 10~^, which could be further improved to sin^ 2^13 < 5.26 x 10~^ in the 
event that ten times larger statistics becomes feasible. 

6 Discussions and Conclusions 

In previous papers [22, 23] we have shown that a Beta-beam experiment where and Pg beams 
produced Tising and ^Li ions respectively are sent from CERN to the magnetized iron calorimeter 
ICAL at INO results in tremendous sensitivity to ^13 and the mass hierarchy. This unprecedented 
sensitivity stems from a combined effect of the CERN-INO distance being magic, as well as the 
energy of the Beta-beam being conducive to encountering near-maximal matter effects, thereby 
resulting in substantial enhancement of the oscillation probability and hence statistics. However, 
while the baseline being magic and the beam energy being in the multi-GcV regime is a virtue for 
probing ^13 and the mass hierarchy, it becomes a serious drawback for CP violation studies. This 
necessitates the deployment of a second baseline option with a lower energy neutrino beam for 
investigating all the three issues simultaneously. The optimal baseline and ion source option for 
CP violation studies was shown in [41] to be L = 600 — 900 km and the ^^Ne and ^He combination 
with intermediate values of 7. In this paper we studied the sensitivity reach of the combined data 
sets from the CERN to INO magical Beta-beam set-up and a CERN to Gran Sasso Beta-beam 
experiment with ^^Ne and ^He as source ions. We reiterate that even for the Neutrino Factory 
experiment, at least two baselines are needed in order to optimally address all the three neutrino 
oscillation parameters. 

For the CERN to INO set-up, the basehne is L = 7152 km and we used ^B and ^Li as sources, 
with boost 7 = 650. The far detector was taken as ICAL, a 50 kton magnetized iron calorimeter 
with detection efficiency of 80%, charge identification efficency of 95%, energy threshold of 1 GcV, 
energy resolution of 0.15£' and neutral current background fraction of 10^^. We call this set-up 
CERN-ICAL@INO. For the CERN to Gran Sasso sector, the basehne is L = 730 km and we 
chose ^®Ne and ^He as source ions. The Lorentz boosts for the ^^Ne and ^He ions are taken 
to be the same and allowed to vary between 250 and 650. Because of the lower beam energy 
for this case it is preferable to opt for a Totally Active Scintillator Detector here. We assume 
an active detector mass of 50 kton with 80% detection efficiency for muons and 20% detection 
efficiency for electrons, energy threshold of 0.5 GeV, energy resolution of 0.03^i?(GeV) for muons 

and 0.06-^ii^(GeV) for electrons and neutral current background fraction of 10^^. We call this 
set-up CERN-TASD@LNGS. We show results for "standard luminosity" where we use the above- 
mentioned detector sizes and efficiencies and 5 x (1.1 x 10^®) and 5 x (2.9 x 10^®) useful ion decays 
in neutrino and antineutrino modes respectively. We also show results for a situation where the 
statistics is ten times larger. 

We probed the physics potential of this two-detector set-up with a Beta-beam as the neu- 
trino source. We have presented our results in terms of three performance indicators. The 
essential results are summarized in Table 4. For 5cy> sensitivity we showed the sin^2^i3 (true) 
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Set-up 


Mass Ordering (3cr) 

NH (True) 


CP Sensitivity (3<t) 
NH (True) 


sin^ 2^13 Sensitivity (30-) 


1.1 X W 

& 

2.9 X 10^^ 


1.1 X W 

& 

2.9 X 10^9 


1.1 X W 

& 

2.9 X 10^^ 


1.1 X W 

& 

2.9 X 10^9 


1.1 X 10^« 
& 

2.9 X 10^^ 


1.1 X 

& 

2.9 X lO^'' 


CERN-INO 
7 = 650, 7152 Km 


4. / X lU 

(4.9 X lO-'') 


n /I w 1 n— t) 

y.4 X lU 
(1.2 X 10-4) 


iNOL 

possible 


iNOt 

possible 


1.14 X 10-3 


1.76 X 10-4 


CERN-LNGS 
7 = 575, 730 Km 


3.89 X 10"^ 
(9.23 X 10-3) 


1.58 X 10-3 
(4.48 X 10-3) 


1.6 X 10-4 
(1.8 X 10-4) 


1.97 X 10-5 
(2.03 X 10-5) 


1.78 X 10-3 


8.59 X 10-5 


CERN-BOULBY 
7 = 575, 1050 Km 


2.49 X 10-3 
(7.87 X 10-3) 


2.19 X 10-4 
(4.1 X 10-3) 


1.85 X 10-4 
(2.02 X 10-4) 


1.99 X 10-5 
(2.04 X 10-5) 


1.41 X 10-3 


1.45 X 10-4 


CERN-LNGS 
7 = 575, 730 Km 
+ 

CERN-INO 

7 = 650, 7152 Km 


2.7 X 10-'' 
(3.58 X lO-'') 


4.64 X 10-5 
(5.45 X 10-5) 


1.42 X 10-4 
(1.49 X 10-4) 


1.78 X 10-5 
(1.88 X 10-5) 


5.46 X 10-4 


5.26 X 10-5 


CERN-BOULBY 
7 575, 1050 Km 

+ 

CERN-INO 
7 = 650, 7152 Km 


2.67 X 10-4 
(3.37 X 10-4) 


4.57 X 10-5 
(5.17 X 10-5) 


1.63 X 10-4 
(1.76 X 10-4) 


1.8 X 10-5 
(1.87 X 10-5) 


6.1 X 10-4 


6.69 X 10-5 


Optimized Neutrino Factory set-up with two improved golden detectors (50 kton each) placed at 
L = 4000 km & 7500 km respectively. = 20 GeV & total 5 x 10^^ decays for ^u" & /i+ each. 


Optimized 

Ncnitriuo F;u-tor\- 


4.5 X 10-5 

(lOO'X ot'r^cp (truo) cov-orai^x") 


1.5 X 10-5 


4.5 X 10-5 



Table 4: Performances of various Beta-beam set-ups at 3(7 in addressing the key unsolved issues: 
mass ordering, CP violation and sin^ 2(^i3 sensitivity reach. For CP sensitivity and mass ordering, 
the minimum values of sin^ 26i^ (true) required for a positive conclusion are presented. Results are 
shown for a five-year run with the reference luminosity: 1.1 x 10^^ (2.9 x 10^^) useful ion decays per 
year in the v [u) mode as well as one order of magnitude higher statistics. The numbers without 
(with) parantheses correspond to 5cp (true) = 90° (5cp (true) = 270°). Note that the sin^ 2^13 
sensitivity reach is independent of the value of 5cp (true) and the true mass ordering because the 
prospective "data" have been generated at 6'i3 = 0. The CERN-INO baseline is insensitive to 5cp- 
For comparison, the expectations from an optimized two-baseline Neutrino Factory set-up with 
upgraded magnetized iron detectors are also listed [4, 40]. 
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reach of the experiment to distinguish maximal CP violation from a CP-conserving scenario. 
We displayed the results for the CERN-TASD@LNGS experiment by itself and when data from 
here are combined with that from CERN-ICAL@INO. From CERN-TASD@LNGS alone we ob- 
tain magnificent sensitivity to Sqp- For NH true and 7 = 575 for ^^Ne and ^He and with the 
standard luminosity, maximal CP violation can be established at 3a at CERN-TASD@LNGS if 
sin^ 26*13 (true) > 1.6 x 10~^ for 6cp (true) = 90°. With ten times more statistics this limit gets 
pushed to sin^ 2^13 (true) > 1.97 xlO^^. Addition of the CERN-ICAL@INO data does not improve 
these hmits significantly. We also explored the dependence of maximal CP violation sensitivity of 
CERN-TASD@LNGS to the detector characteristics of TASD. Resuhs are similar if IH is true. 

Measurement of ^^13 and mass hierarchy can be performed extremely well at CERN-1CAL@1N0. 
We studied the impact of adding the data from CERN-TASD@LNGS on the final combined 
sensitivity to these parameters. For the mixing angle ^^13 we defined the sin^ 2^13 sensitivity as 
the range of sin^ 2^13 which could fit at 3a the data generated for sin^ 2^13 (true) = 0, after 
marginalization over all oscillation parameters and the density profile. We found that the sin^ 2^13 
sensitivity of both set-ups taken alone are comparable. CERN-ICALOINO could limit sin^ 2^13 < 
1.14 X 10"^ for the standard luminosity and with 7 = 650, while CERN-TASD@LNGS could 
restrict sin^2^i3 < 1.78 x 10~^ with 7 = 575. For the enhanced luminosity these hmits would be 
sin^2^i3 < 1.76 X 10~^ and sin^ 2^13 < 8.59 x 10~^ respectively. However, when added together, 
the combined data can limit sin^ 2^13 < 5.46 x 10~^ for the standard luminosity and sin^ 2^13 < 
5.26 x 10^^ with a factor of ten enhanced luminosity. 

The sensitivity to sgn{Am'^i) is defined in terms of the sin^ 26*13 (true) required in order to rule 
out the wrong hierarchy at 3a. We had shown earher that for CERN-ICAL@INO with 7 = 650 
and standard luminosity the wrong inverted hierarchy can be disfavored at 3a if sin^ 2^13 (true) > 
4.7 X 10-"^ for Sep (true) = 90°. With ten times more statistics this improves to sin 2^13 (true) > 
9.4 X 10"^ We found in this paper that the sgn{Aml^) sensitivity of CERN-TASD@LNGS alone 
is rather poor in comparison. However, when we add the two data sets, we find an enhancement in 
the sensitivity. With data from both set-ups taken together the wrong inverted hierarchy can be 
disfavored at 3a if sin^ 2^13 (true) > 2.7x 10~^ with the standard luminosity when Sqp (true) = 90°. 
With ten times more statistics it would become possible if sin^ 2^13 (true) > 4.64 x 10~^. Results 
for Scp (true) = 270° are similar, and so arc the sensitivity reaches for IH true. 

Another choice for an intermediate baseline could be sending a beam from CERN to the Boulby 
mine in UK. The CERN to Boulby distance is about 1050 km. Since INO and Boulby mine are in 
opposite hemispheres with respect to CERN, it might be easier to contemplate a decay ring design 
which could be used to send beams to both these (plausible) detector sites. While in the text we 
explicitly presented and discussed results for the CERN-LNGS baseline, Table 4 shows results for 
the CERN-Boulby mine set-up as well. All detector characteristics are taken to be the same as 
for the LNGS case. The sensitivity reach can be quantitatively seen to be comparable for both 
the intermediate baseline options when combined with CERN-INO data, with the CERN-LNGS 
baseline being marginally better. 

The sensitivity obtained in the set-ups considered here could be compared to that obtained 
in a benchmark high 7 set-up considered in [33]. In this paper, the authors proposed using a 
7 = 350 ^^Ne and ^He Beta-beam sent to a megaton water detector located at a distance of 730 
km from CERN. For total exposure of 5.0 Mton-year, this set-up returns a 3a sin^ 2^13 sensitivity 
of 5.7 X 10"^. Normal mass hierarchy discovery reach at 3a ranges between 2.4 x 10~^ and 
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1.6 X 10 ^ depending on the choice of 6cp (true). Maximal CP violation can be established at 
3(j if sin 2^13 (true) > 5.2 x 10"^ for Scp (true) = 7r/2 and sin2 2^i3 (true) > 5.5 x 10"^ for 
5cp (true) — 37r/2. 

The Beta-beam set-ups discussed here have a reach comparable to those obtained from exper- 
iments based on optimised Neutrino Factories, see Table 4. Here, the reference Neutrino Factory 
set-up we have considered is the version with two "golden magnetized detectors" of 50 kton each 
at L = 4000 km and 7500 km, a detector threshold of 1 GeV, 5 years n^/n^ data with — 
20 GeV and a luminosity of 10^^ decays per year for each polarity, with a background fraction of 
5 X 10~^. Further details of these set-ups and their resolutions and efficiencies can be found in 
[4, 40]. 

It is of interest to compare the sensitivity reach of our two-detector Beta-beam set-up with 
that chosen in [31], where the authors used L — 2000 km for the shorter basehne magnetized iron 
calorimeter detector in addition to another at the magic baseline, and ^B and ^Li ions as sources 
for both. While they use a smaller boost factor of 350 for the beam, they require larger number of 
ion decays than we have used for our standard luminosity. Our results with standard luminosity 
are comparable to those of [31] with larger luminosity, as noted in the Introduction. Similar reach 
for CP sensitivity in the two proposals even though with different source ions and intermediate 
baseline is not surprising. The relative optimization between 7 and luminosity for the two sets 
of sources was studied in detail in [41]. It was shown that to obtain the same "physics" with 
^B/^Li and ^^Ne/^He, one needs 12 times larger luminosity for the former and 3.5 times larger 
7 for the latter. While we have used a factor of 1.86 higher boost than that in [31], the authors 
of [31] consider 7.39 (2.80) times higher luminosity for their neutrino (antineutrino) mode. For 
lower number of ion decays the sensitivity of the set-up proposed in [31] becomes worse than our 
set-up. The main strength of the set-up proposed here emerges from the fact that use has been 
made of the optimal choice for the source ions and baselines [41] for CP, sgn{Am1i) and 6*13. 
The boost factors that we use, though larger, should be plausible at CERN and were optimized 
for the CERN-INO set-up for determining sgn{Aml^) and ^13 [23]. For the corresponding boost 
possible for ^^Ne and ^He we achieve good sensitivity to CP violation in the CERN-LNGS set-up. 
The combined data set from both detectors provide exceptional sensitivity to all three neutrino 
parameters. With ten times larger luminosities, which are also considered to be plausible [43], the 
sensitivity of our set-up escalates to unprecedented levels and arc similar to those possible with a 
high performance two-detector Neutrino Factory set-up, using improved detectors with very low 
backgrounds and low threshold (see Table 4) . 

In conclusion, we show that a Beta-beam experiment with two carefully chosen detectors at 
optimal distances can provide unprecedented sensitivity to the establishment of CP-violation in 
the lepton sector, addressing the issue of neutrino mass ordering, and determination of the mixing 
angle 9 13. We have considered a Totally Active Scintillator Detector at the shorter baseline of 730 
km and a magnetized iron calorimeter at a near-magic baseline distance of 7152 km, both of mass 
50 kton. In tandem, they provide an exceptional sensitivity to the neutrino parameters and match 
the precision achievable in high performance Neutrino Factory set-ups. 
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A Appendix: Background Rejection 

In this Appendix we discuss some technical issues pertaining to rejection of the beam related and 
atmospheric neutrino backgrounds. 

A.l Beam related neutral current backgrounds 

Prom Fig. 1, one can see that at the CERN-LNGS baseline (L = 730 Km), the expected event 
rate sharply depends on the choice of Scp and ^13. Assuming 1.1 x 10^^ useful decays per year 
for ^^Ne and at sin^ 29i3 = 10^^ with a boost factor of 575 and normal hierarchy, the variation 
in the signal in 5 years with 5cp has been depicted in Table 5. The beam related neutral current 
background in this case (after folding with 10~^ suppression factor for TASD) comes out to be 23 
events. So in this case, the background over signal ratio may vary over a wide range from 0.26 to 
1.92. This case is quite different to that encountered in the CERN-INO Beta-beam experiment 
[22, 23] where we hardly expect any beam related background. Though the total signal to neutral 
current background ratio is high, we stress that an aspect which is extremely crucial is the different 
spectral shape of the signal compared to the background. In our analysis we have cautiously taken 
into account the beam related background as described in the appendix of [23]. 



^CP 





90° 


180° 


270° 


Signal (5 yrs) 


67 


90 


35 


12 



Table 5: Number of u signal events in 5 years for various 5cp with a 50 kton TASD detector 
(sin^ 2^13 = 10-3 and NH with 7 = 575). 

A. 2 Atmospheric backgrounds 

Next we turn to the backgrounds due to atmospheric neutrinos. In the CERN-LNGS set-up, we 
have considered a TASD detector with an energy threshold of 500 MeV. With ^^Ne and 7 = 575, 
we have an energy spectrum extending upto 3.9 GeV with a peak at 2.25 GeV. In this energy 
range the number of muons from atmospheric neutrinos giving tracks inside the detector along 
the direction of the /3-beam flux are expected to be very large (see Table 6). Without demur one 
can see that the atmospheric backgrounds are huge compared to the signal rates. Also the TASD 
detector does not have any charge identification (CID) capability and it makes the case more 
arduous in the sense that we have to consider both /x" and /i"*" atmospheric events simultaneously 
as the backgrounds when we expect a signal from the neutrino beam alone. In such a situation, 
the timing information of the ion bunches inside the storage ring turns out to be a valuable tool 
to tackle these backgrounds. 
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Energy range (GeV) 


/X events 


events 


0.5 - 4 


1342 


447 


1 - 4 


852 


344 



Table 6: Expected atmospheric neutrino events along the beam direction in 5 years with a 50 kton 
TASD detector. 

A. 3 Bunch size and background reduction 

For a 5T magnetic field and 7 = 575 for ^^Ne ions, the radius of the curved section would be 
R ~ 720 m. The useful decays are those which occur in one of the straight sections where emitted 
neutrinos fly towards the direction of the detector. If is the length of one of the straight sections 
of the decay ring and if we demand that L^/ (27ri? + 2Lr) ~ 36% (useful fraction of ion decays) 
then Lr comes out to be 5818 m. With this design, the total length of the storage ring {27rR + 2Lr) 
would be 16161 m. Now if we can tolerate the effective atmospheric background over beam signal 
ratio to upto 10~^ then 



10 



_^ Atmospheric Background x Suppression Factor (Sf) 

Beam Signal 



(A-i; 



The total atmospheric background in the energy range of 0.5 to 4 GeV is 1789 (see Table 6) and 
the expected lowest signal rate in 5 years from the ^^Ne /9-beam is 12 (with sin^ 2^13 = 10~^ and 
Sqp = 270°). So from Eq. A-1, wc have Sf = 6.7 x 10^^. To achieve this Sf, one has to ensure an 
excellent timing information for the signal correlating with the source bunches. This requirement 
puts a constraint on the maximum allowed time-length of the ion bunch Tf, inside the storage ring. 
The Sf and Ti, are related in the following fashion. 



S, = X 0.36 



^ 2nR + 2Lr 



(A-2) 



where v ~ c is the velocity of the ion and A*";, is the maximum number of bunches of the ion 
circulating inside the storage ring at the same time. One can readily see from Eq. A-2 that 
Sf — 6.7 X 10~^ can be achieved with 12.5 ns time- length of the ion bunch taking A^^ = 8. 

If we increase the threshold energy of the TASD detector from 0.5 to 1 GeV then we can 
suppress the atmospheric events by a factor of 1.5 (see Table 6) and now the total atmospheric 
background would be 1196. But in this process, the signal rate also gets reduced from 12 to 8, 
keeping the Sf almost same. 

Another important point one should keep in mind is that if we can work with one order higher 
luminosity of the Beta-beam flux then the signal rate will be enhanced by a factor of 10 and 
therefore one can relax the time length of the ion bunch from 12.5 ns to 125 ns provided that Nt, 
is kept fixed. 
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